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Abstract 

In this note, I present a simple approximate expression for the number C(3), 
known as the Apery's constant, which is accurate to 21 digits. This finite 
closed-form expression has been found experimentally via the PSLQ algo- 
rithm, with a suitable search basis involving the numbers vr. In 2 , In (1 + -\/2 ), 
and G (the Catalan's constant). The short Maple code written for finding 
the rational coefficients is also shown. 
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"/n fact, numerical experimentation is crucial to Number Theory, 
perhaps more so than to other areas of mathematics. " 

F. R. Villegas 

1. Introduction 

The Apery's constant is defined as the number to which the series Xl^i ^/^^ 
converges and it is so designated in honor to R. Apery, who proved in 1978 



that this number is irrational It is not known if it is a transcendental 
number. The convergence of that series to a real number between 1 and 2 
is guaranteed by the Cauchy's integral test, a result that remains valid for 
l/'T-*, for any real s > 1. In fact, convergence is also found for > 1, 
a complex domain in which this series is defined as C{s)i the Riemann's zeta 
function. Apery's constant is then identified with (^(3) = 1.202 056 903 . . . 

For positive integer values of s (s > 1), it was Euler (1734) the first to 
obtain an exact closed-form result for C,{s)i namely C(2) = 7r^/6, the solution 
of the so-called Basel's problem |2|. Some years later, he did generalize this 
result for all (positive) even values of s [s], 0]: 

o2n— 1 r) ^2n 

where n is a positive integer and G Q are Bernoulli numbers. For odd 
values of s, closed-form expressions are not known. In a recent work, it is 
even claimed that C{2n + 1) is not a rational multiple of vr^"'"'"^ 23|. 

The increase of interest in C(3), which comes from both pure and applied 
mathematicsjll has stimulated its high-precision numerical computation, as 



^For instance, given three integers chosen at random, the probabihty that no common 
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well as the search for simple approximate expressions for it jsl- Let us adopt 
a reasonable criterium for the adjective "simple," in the context of finite 
approximate expressions. Here, this will designate closed-form expressions 
containing only a few terms/factors composed by other mathematical con- 
stants and a few integers with small absolute values. This criterium is, of 
course, vague due to the forms "a few" and "small". To make it not-so- 
imprecise, "a few" will mean less than, say, ten (certainly not twenty), and 
"small" will mean less than, say, five hundreds (certainly not a thousand). 
A such approximation was presented by Galliani (2002), namely 6| 

-3 



C(3)^^ + ^ (1 + 27- , ) , (1) 

^ vr* \ 130 + TT^ ' 



where 7 is the Euler-Mascheroni constant, which is accurate to 4 digits (typo 
corrected). Another nice simple approximation is 



C(3)«^7+§. (2) 

due to Hudson (2004), which is accurate to 7 digits j6|. Among the many 
approximations for C(3) presented by Hudson, the most accurate is ^ 

C(3) ^ 525587 1/^^, (3) 

which is accurate to 12 digits and, clearly, is not a simple approximate ex- 
pression (in our terminology). The same for C(3) ~ ^H^^^ vr^, which I found 



factor will divide them all is 1/C(3). Also, if n is a power of 2, then the number 
of distinct solutions for n = p + xy with p prime and x, y positive integers obeys the 
asymptotic relation ^{n)/n ^ 105 C(3)/ (Stt"*). It also arises in a number of physical 
problems, including the computation of electron's gyromagnctic ratio. 
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by searching for an integer relation between ({3) and vr^, before the appear- 



ance of Ref. 23|. In trying to reach greater accuracy, however, one soon 
observes that it is very difficult to avoid the appearance of large integers. 

In this short note, I present a simple approximate expression for ({3) 
which is accurate to 21 digits, thus triply most accurate then the best pre- 
vious simple approximation. This closed-form expression involves only small 
integers and the constants vr, ln2 , In (1 + a/2), and G (the Catalan's con- 
stant) and it has been found experimentally by employing a search basis with 
six elements in the PSLQ algorithm. The short Maple code I have written for 
finding such expression is also shown in view to stimulate people to develop 
their own computational experiments. 

2. The search for integer relations and the PSLQ algorithm 

An important task in experimental mathematics is to search for inte- 
ger relations involving a finite set of computed numbers. An integer re- 
lation algorithm is a computational scheme that, for a given real vector 
X = {xi,X2, . . . , Xn), n being a positive integer, n > 1, it either finds a nonnuU 
vector of integers a = (ai, 02, . . . , a„,) such that ai a;i + 02 X2 + . . . + a„, x„ = 
or else establishes that there is no such integer vector within a ball of some 
radius about the origin. 

Presently, the best algorithm for detecting integer relations is the PSLQ 
algorithm (acronym for Partial Sum of Least sQuares) introduced by Fer- 
guson and Bailey (1992) A simplified formulation for this algorithm, 
mathematically equivalent to the original one, was subsequently developed 



by Ferguson and co-workers (1999) 



This more efficient version of PSLQ 
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is currently implemented in both Maple and Mathematica, two of the most 
popular mathematical softwares. It was this version, together with some 
reduction schemes, that was named one of the 'ten algorithms of the cen- 



tury' IqI. 

In short, PSLQ operates as follows. Given a vector x of n given real 
numbers, input as a list of floating-point (FP) numbers, the algorithm uses 
QR decomposition in order to construct a series of matrices Am such that the 
absolute values of the entries of the vector = ■ x decrease monoton- 
ically. At any given iteration, the largest and smallest entries of ym. usually 
differ by no more than a few orders of magnitude. When the desired in- 
teger relation is detected, the smallest entry of y.^ abruptly decreases to 
roughly the computer working precision e and the relation is given by the 
corresponding column of A^. This numerically stable matrix reduction pro- 
cedure, together with some techniques that allow machine arithmetic to be 
used in many intermediate steps, usually yields a rapid convergence, which 



10| . When a relation 



makes PSLQ faster than other concurrent algorithms 
is detected, the ratio between the smallest and the largest entry of the vec- 
tor A~^ ■ X can be taken as a "confidence level" that the relation is true 
(i.e., exact) and not an artifact of insufficient precision. Very small ratios at 
detection certainly indicate the result is probably truej^ 



^In addition to possessing good numerical stability, PSLQ is guaranteed to find an 
integer relation in a number of iterations bounded by a polynomial in n. 
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3. Maple code for finding an approximate expression via PSQL 

Since an efficient PSLQ routine is available as part of a Maple package 
named IntegerRelations, then simple short codes can be written in this 
system. For illustrating this, let me list just the source code I have written 
for finding an approximate expression for C(3). 

> restart; # Clear memory 

> Digits := 24: # The number of digits for FP numbers 

> with(IntegerRelations) : # Call the package containing PSLQ 

> xSymb: = [Zeta(3) , 1 ,Pi''2*ln(2) ,Pi*ln(2) "2, ln(2) "3,ln(l+sqrt (2) ) "3,Pi*Catalan] ; 

xSymb := [C(3),l, 7r2ln(2), 7rln(2)2, ln(2)3, ln(l + ^/2f, vr Catalan] 

> n := nops(xSymb); # The number of elements in xSymb 

n := 7 

> X : =evalf (xSymb) : # Convert to FP numbers 

> a:=PSLQ(x); # Applies PSLQ algorithm to x 

a := [10, 394, -11, 283, -472, -209, -186] 

> soma:=0: 

> for k from 1 to n do soma:=soma+a[k] *xSymb[k] ; end do: 

> aprox : =solve (soma=0 , Zeta(3)); 

5 11 , , 283 , , 236, , ,o 209, , 93 tt Catalan 

aprox := \ tt^ ln(2) 7rln(2)2 H ln(2f H ln(l + V2f H 

^ 197 394 ^ ^ 394 ^ ^ 197 ^ ^ 394 ^ ^ 197 

> evalf ( aprox, 30 ) ; # Our approximation 

1.20205690315959428539958993430 

> evalf ( Zeta(3) , 30 ) ; # Exact value (for a better comparison) 

1.20205690315959428539973816151 
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With this Maple routine, I have found the following simple approximate 
expression for the Apery's constant: 

5 11 2 , o 283 , 2o 236 , . 209 , 3 93 ^ 

((3) ^ \ TT^ In 2 TT In^ 2H In^ 2H In^ (1 + V2 )H vr G, 

^ 197 394 394 197 394 ^ ^197 

(4) 

which is accurate to 21 digits, as the reader can check by comparing the last 
two outputs of the program. Equation (jlj) is the main result of this paper. 
All that rests is to explain the motivation for choosing the vector 

X = (^1, vr^ ln2, vr \n^2, ln^2, In^ (1 + ^2), n 

as the search basis for C(3)- The motivation comes primarily from a conjec- 
ture by Euler (1785) that C(3) = a In^ 2 + /3 ttMu 2 , for some a, /3 G Q (l2|. 



After many numerical experiments, I could not find any such pair of rational 
coefficients for composing an exact closed-form expression for C(3)- Not even 
a simple approximate expression was found with only these two terms. I am 
currently testing the conjecture by Connon (2008) that either a or could 
contain a factor of \/2, or another small surd, or perhaps In (2-7r) 13|. I was 
then inclined to improve the basis by including other weight-3 constants such 
as vr In^ 2 and ttG, which arise in many exact results involving C(3), as for 
instance y 



(i) Special values of non-elementary functions: 



a3)-^7r2 In2+|ln32, [Lis (if^)] = f C(3) + ^ ln^2- 
isl, 7r2V;(-^)(l) = Ivr^lnA+^vr^ In2+^7r2 ln7r+iC(3) U 



Lis (I) = 
ife^' ln2 
(typo corrected); 

' Weight, here, is according to the definition by Bores and Moll (see Ref. jl^, p. 203). 



(ii) Series: 



E 



oo 



(-1) 



fe+1 



f C(3)-47rG 



iC(3) 

'^'^ Em=0 



iln^2[l^j, 



^ Ea 



iTT^-ivr^ ln2- 



oo C(2fc) 
fc=l (fe+1) 24*= 

oo (-ir^^-Wi(i) ^ i^2i^2-|C(3) [19| 



(2m+5)! 



(iii) Definite integrals: 
J^^"^ X In (sin x) dx 



IttG-^tt^ ln2 



^C(3)-| 



In 2 , JJ^ X In (cos x) dx = — ^ C(3) + 
20|, /g'^^x^ tanxcix = -^C{3) + \tt G-j^rr^ In 2 y, 



/o"/^ xVsinxrfx = 27rG'-|C(3) 
(see Eq.(6.6.25) in Ref. Q). 



/c 



1 (arcsint)-^ 







dt=i7r2ln2-|C(3) 



Here, Li3(2;) := E^^i z'^/n^ is the trilogarithm function, iIj^~^\z) is tlie poly- 
gamma function (extended to negative indexes, according to Ref. (2l|), and 
E2m+i{x) are Euler polynomials. 

With respect to In^ (1 + V2), which is the cube of a logarithm of an 
algebraic number, thus a weight-3 constant, I have included it in the search 
basis because, curiously, the number arcsinh(l) = In (1 + a/2) emerges in 
the coordinates of the vertices of a cusped hyperbolic cube whose volume is 
I C(3), as I have found in a recent work on multiple integrals related to C(2) 
and C(3) [221 . 
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